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Abstract 

We investigate the form of the non-perturbative parameterization in both the impact 
parameter (b) space and transverse momentum (jot) space resummation formalisms for the 
transverse momentum distribution of single massive bosons produced at hadron colliders. We 
propose to analyse data on T hadroproduction as a means of studying the non-perturbative 
contribution in processes with two gluons in the initial state. We also discuss the theoretical 
errors on the resummed Higgs transverse momentum distribution at the LHC arising from 
the non-perturbative contribution. 



* Anna.Kulesza@bnl.gov 
'W.J. Stirling® dur ham. ac .uk 



1 



1. Introduction 



The main discovery channel for a light Higgs boson at the future Large Hadron Collider (LHC) 
is the gluon-gluon fusion process gg — ► HX. Studies of the production characteristics are of 
great importance for further understanding QCD and determining the Higgs parameters. A 
particularly interesting quantity is the transverse momentum (pr) distribution of the produced 
Higgs boson. In particular, a precise knowledge of the Higgs pr distribution is important for 
determining an optimal set of cuts on the final state particles. In the past, the study of the 
transverse momentum distribution in Drell-Yan lepton pair and electroweak boson production 
helped to establish the validity of the soft gluon resummation formalism for the distribution at 
small pt, and played an important part in determining the Standard Model (SM) electroweak 
parameters. 

Reliable predictions for transverse momentum distributions at small px can only be obtained 
if soft gluon emission effects are correctly taken into account. The soft radiation manifests itself 
in the presence of large logarithmic corrections in the theoretical expressions. The two prominent 
examples are recoil logarithms of the form \ a i 2n ^ 1 ) (Q 2 /p' T ) and threshold logarithms of the 
form of In*' 2 ™" ^ ( 1 — z ) I ( 1 — z) (where z = Q 2 /s). A finite result can be recovered by resumming 
these corrections to all orders in as- To date, the most studied method for resumming the recoil 
(or 'Sudakov') corrections in transverse momentum distribution for vector boson (7*, Z, W) 
production is the impact parameter (b, Fourier conjugate to pr, space) formalism, also known 
as the Collins-Soper-Sterman (CSS) formalism [1]. Derived from the b space resummation 
formalism are pt space methods [3, 4, 5], which not only provide a very good approximation 
of the b space result but also avoid certain drawbacks related to the b space method [3, 4]. 
As recently demonstrated, both the recoil and threshold corrections can be taken into account 
within the framework of the joint formalism [6, 7, 8]. 

The application of the CSS formalism to Higgs production at the LHC is a straightforward 
task and several analyses exist in the literature [9, 10, 11, 12]. One of the topics discussed here 
is the application of the pt space formalism to Higgs production at the LHC. Some preliminary 
results based on the pt space method have already been presented in Ref. [13]. 

The CSS formalism resums logarithms arising from arbitrarily soft gluon emission and there- 
fore needs a prescription for dealing with the Landau pole. Originally this was incorporated by 
introducing an arbitrary (inverse energy) scale (&*) above which the perturbatively calculated 
expression remained constant, i.e. 'frozen' at 6*. Additionally, a gaussian function, for example 
exp(— gb 2 ) in its simplest form, was introduced to correct the formalism for non-perturbative 
effects at large b, i.e. intrinsic transverse momentum. The parameters of this function are in 
principle determined from experiment through fits to the data. In the case of Drell-Yan lepton 
pair or electroweak boson production, determining the size of the non-perturbative contribution 
associated with the incoming quarks can be done relatively precisely since there is a significant 
amount of experimental data available. 

Other methods for dealing with the Landau singularity have also been proposed in the 
literature. The approach of [6, 7], initially developed for the joint formalism and later applied 
to the b space formalism in [12], relies on performing the inverse Fourier transform from b to pr 
space as an integral along a contour in the complex b plane. It returns a well-defined resummed 
distribution for all nonzero values of px, even without introducing any extra non-perturbative 
function. It turns out, however, that some non-perturbative input is still needed to obtain a 
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good description of the Z production data at the Tevatron [7] . In another development of the 
b space formalism, Qiu and Zhang [27] proposed a smooth extrapolation of the perturbative 
resummed cross section to push the Landau singularity to infinity. 

A complete description of Higgs production within a resummation formalism also requires 
a prescription for dealing with the non-perturbative regime. While the 6* prescription (or 
other prescriptions) for defining the non-perturbative regime can be kept the same, a different 
non-perturbative function is expected to describe intrinsic kx in the gg channel. A common 
hypothesis is that the non-perturbative parameters for this (gg) channel can be determined by 
multiplying one of the parameters of the non-perturbative function for Drell-Yan (qq) production 
by a factor of Ca/Cf, although to the best of our knowledge there is no rigorous theoretical 
basis for this procedure. 

In this paper we calculate the Higgs transverse momentum distribution at the LHC using 
the pr space formalism. This method also needs a non-perturbative (gaussian) input function to 
account for the intrinsic kx- In order to determine the non-perturbative coefficients we propose 
here to study transverse momentum distribution in T production process. We argue that the 
(gg — > T) hadroproduction provides valuable information on the amount of intrinsic kx carried 
by incoming gluons. The analysis begins by testing the existing forms of the non-perturbative 
parameterization for Drell-Yan production and the Ca/Cf hypothesis in b space. A similar 
study is then performed for the px space formalism. We finish by considering (or, in the case of 
b space analysis, re-examining) the effects of the non-perturbative function on Higgs transverse 
momentum distribution at the LHC, and draw some conclusions on the uncertainty in the 
prediction for the Higgs boson transverse momentum distribution arising from non-perturbative 
effects. 



2. Theory 

We begin by recalling basic formulae for both b space and px space formalisms. A brief review 
of non-perturbative parameterizations follows. 

2.1. b space 

The b space resummed part of the theoretical cross section for Higgs boson production through 
the gluon-gluon fusion in hadronic collisions p p — > H + X has the following form, cf. [1] 

dp 2 dQ 2 = °~° ^ ~ m H) y o dx A dx B o \ xaxb —I x 

^jf dbb J (p T b) exp[S(h,Q)]F NP (b 1 Q,XA,x B )f' g/ A ^) f'g/B (^b, , 

with a = ^ G 5 jq^ — , t = Q 2 /s, b = 2exp(-7 £ ), and where 
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For the gluon-gluon fusion process [14, 15], 
= 2C A , = 
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and fa = fC A - §N F T R . 

The 'modified' parton distributions /' are related to the MS parton distributions, /, by a 
convolution [1, 16, 17] 



fg/H( x >») =Y1 j x ^ C 9 a {~^) fa/H(z,V) , 
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The next-to-leading logarithm (NLL) accuracy requires knowledge of the , B^\A^ and 
Cg? coefficients; the coefficients A^ , i?( 2 ) and Cga are of NNLL order. 

The expression (1) provides a good description of the pr distribution at small pt] for larger 
values of p? one needs to match the resummed result with the fixed order result. This is achieved 
by writing 

da da TCS 

where Y(pt, Q) is the difference between the fixed-order and resummed results expanded up to 
the order at which the fixed-order expression is considered. 
The 6* variable in (1) is defined as 

b 



V 7 ! + (W ' 



6* < 



(7) 



which ensures that the resummed b space expression is well-defined. The replacement of b by 6* 
in (1) prevents the argument of as from entering the non-perturbative regime by 'freezing' the 
perturbative contribution at a certain b ~ b\ im . 

Using renormalization group analysis arguments, in [1] a universal form of the non-perturbative 
function F/vp was proposed: 



F NP (Q,b,x A ,x B ) = exp 



-h Q (b)\n 



_Q_ 

2Q C 



hA(b,XA) ~ h B (b,XB) 



(8) 



where Qq is an arbitrary constant indicating the smallest scale at which perturbation theory is 
reliable, Qo ~ l/&ii m - The functions Iiq, h a , h\, are to be extracted by comparing theoretical 
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predictions with experimental data. As postulated in [1], the flavour dependence of F NP can 
be ignored. The ln(Q/Qo) dependence in (8) is required to balance the Q dependence of the 
Sudakov factor exp 5. In addition, hn was proved to be universal and its leading b 2 behaviour 
at large b is suggested by the analysis of the infrared renormalon contribution [21], as well as a 
recent analysis of the dispersive approach to power corrections [22]. 

However, the detailed form of the non-perturbative function F^ P (Q, b, x a ,Xb) has remained 
a matter of theoretical dispute. Early studies of fixed-target Drell-Yan experimental data, see 
e.g. [19], suggested that a Gaussian parameterization of an intrinsic qx distribution provided a 
good description of data in the low pt (1 — 2 GeV) regime. Motivated by this result, Davies et 
al. (DSW) [20] approximated the function F NP by 



F NP (Q, b, x A ,x B ) = exp 



-^ 2 ln(^-)-^ 



(9) 



The gi parameter in (9) can be interpreted as a measure of the intrinsic transverse momentum, 
whereas 52 represents a contribution coming from unresolved gluons with kx < Qo as the 
structure functions evolve from scales 0(Qo) to 0(Q). Assuming (9), with a particular choice 
of Qo = 2 GeV, b\- im = 0.5 GeV^ 1 , and using the Duke-Owens parton distribution functions [23], 
the DSW analysis gave 

51 = 0.15 GeV 2 , g 2 = 0.40 GeV 2 . (10) 

An alternative parameterization, proposed by Ladinsky and Yuan (LY) [24], incorporates an 
additional dependence on r = x a Xb 



F NP (Q,b,x A ,x B ) = exp 



-g 2 b 2 In (^^j ~ 9ib 2 - gig3bln(l00x A x B ) 



(11) 



Choosing Qo = 1.6 GeV, bu m = 0.5 GeV 1 and using the CTEQ2M parton distribution func- 
tions, the parameters in (11) were determined, 

51 =0.1ll°;°3 GeV 2 , g 2 = 0.58^ GeV 2 , 53 = -1.5±g;i GeV" 1 . (12) 

Both parameterizations were reviewed in [25]. Using high-statistics samples of Drell-Yan and 
(CDF Tevatron Run 0) Z production data, the values of the DSW parameters were updated, 

5i = 0.24t°;°7 GeV 2 , g 2 = 0.34lg;{g GeV 2 , (13) 

and the LY parameters were found to be 

51 = 0.15^^3 GeV 2 , 52 = 0.48±2;°* GeV 2 , 53 = -O.SSlft* GeV" 1 . (14) 

Recently, a new form of Gaussian parametrization with an x-dependent term, proportional 
to b 2 (as opposed to the term linear in b in Eq. 11) has been proposed [26]: 

F NP (Q, b, x A , x B ) = exp \-g 2 b 2 In ( - 9l b 2 + 5153b 2 ln(100x A x B )) 



-.v.j - —j (I5) 

to provide a good description of the low Q Drell-Yan data and the CDF, DO Z data from 
Tevatron Run and Run I. The values of the coefficients determined in [26] are: 

51 = 0.21 ± 0.01 GeV 2 , 52 = 0.68 ± 0.02 GeV 2 , 53 = -0.60±g;g| GeV 2 . (16) 

In our analysis described in the following sections we choose to use the standard 6* pre- 
scription for the CSS formalism. An essentially similar analysis of non-perturbative effects can 
be performed for the CSS formalism with different prescriptions presented in [6, 7, 12] and 
in [27, 11]. 
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2.2. p T space 

The resummed expression in px space, corresponding to (1), has the following form [5] 



da 



dpxdQ 2 



aoT7r5(Q 2 — Mff) J dxA dxs S {^xaxb 
i 1 dpr * , 
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(18) 



and Ei = — px dT<2 /dpr- The factor S v and the c and r coefficients are listed in [4]. 

We choose to incorporate the low energy effects using the form of the pt space non-perturbative 
function F NP advocated in [3], 



■ ? NP 



1 — exp [— dp\] 



(19) 



The role of this function is to account for the distribution in the very low px region, and here we 
are assuming that the shape there is approximately gaussian. However, in order to combine this 
with the perturbative result, the latter needs to be 'frozen' or 'switched off' at some critical value 
of pt where the coupling as becomes large. A similar freezing is required in the b space approach 
where the coupling is effectively as(l/b). In other words, in a similar fashion to the b space 
method, we require not only (i) a form F NP for the distribution in the non-perturbative region, 
but also (ii) a prescription for moving smoothly from the perturbative to the non-perturbative 
region. One possibility for the latter is the 'freezing' prescription of [3], 



Pt +Pt iim ex P 



which has the property 
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It is important to note that there are two pieces of information contained in this definition: the 
value of the limiting value priim and the abruptness of the transition to this value. The use of a 
gaussian function in the definition (20), compared to say a power law function, implies a rapid 
transition from the perturbative to the non-perturbative region. 

As noted earlier, the considerable amount of data on transverse momentum of colour singlets 
produced in processes with two initial quarks makes it possible to determine the non-perturbative 
function relatively precisely. This is not the case for the gluon initiated processes. In fact the 
only data available are on the transverse momentum distribution for low Q resonance production. 
We argue here that E605 data on T hadroproduction [29] can provide useful information on the 
parameters of a non-perturbative function for processes with gluons in the initial state. At the 
values of xf where the cross section is measured, the qq contribution to T production can be 
neglected leaving only the gg channel [28]. Moreover, we also argue that final state interactions 
do not introduce significant additional effects and therefore can be neglected. Since the final 
state interactions are expected to give rise to power corrections, their effects at values of Q ~ my 
relevant for T production are at the percent level of the total cross section. Consequently, the 
basic mechanisms responsible for T and Higgs production are similar. We will therefore assume 
that a good description of the transverse momentum in T production can be achieved with the 
help of the standard formulae for gluon initiated processes in b or pr space, while allowing for an 
arbitrary overall normalization. We also assume that the major contribution to the measured 
cross section for T production comes from the IS resonance with ?^y(is) = 9-5 GeV, and we 
therefore neglect the contributions from other resonances in our theoretical predictions. 

3. Non-perturbative function in b space 

We begin our analysis by studying the non-perturbative function for the impact parameter 
formalism in b space. We assume the standard value of the bn m parameter, i.e. b\\ m = 0.5 
GeV -1 . The dependence of a form of the non-perturbative function on the by mi parameter 
and correlations between the non-perturbative parameters and bn m are not investigated here; a 
discussion of this issue can be found in [17]. 

Given the r-independent form of the non-perturbative function and the relatively narrow 
mass ranges of the data sets under consideration, the fitting procedure can be greatly simplified 
if one adopts an effective form for the non-perturbative function 

f np = e - 9 b 2 ( 22) 

with g an effective parameter, fitted separately for each of the pt distributions in the different 
mass bins. In this way one can test the dependence of the non-perturbative function on Q, and 
in particular possible departures from the proposed logarithmic dependence. Then, from the 
effective values of the non-perturbative parameters g at different values of Q, the values of the 
non-perturbative parameters g\ and gi and/or g% can be deduced. All fitting is done using the 
MINUIT fitting programme [35] . The resummed predictions fitted to data are accurate up to the 
NNLL level, but without including an available numerical estimate for the coefficient [37], 
the effect of which is known to be very small. The Y term in (6) is taken at leading order. Since 
all the data analysed are at small px, the accuracy of matching does not play a significant role 
here. We use MRST2001 parton distribution functions [36] in all our theoretical predictions. 

As a test of the fitting method using the effective function (22), we attempt to reproduce 
the values of the g±, gi parameters obtained by the BLLY collaboration [25]. We take exactly 



6 



the same set of data as in Ref. [25], i.e. the first two mass bins (7 < Q < 8 GeV and 8 < Q < 9 
GeV) of the E605 data [29] (data with p T < 1.4 GeV), R209 data [30] (p T < 2 GeV) and CDF 
Run data on Z boson production [31] (data with px < 23 GeV). The p? distribution measured 
by the R209 experiment in the 8 < Q < 11 GeV bin was read off the plot in [25] . 1 The authors 
of [17] pointed out the impact of the uncertainty in the normalization of the px distributions 
on the determination of the non-perturbative parameters. Here we deal with the normalization 
uncertainty by simultaneously fitting multiple non-perturbative effective coefficients g and one 
common normalization factor to data gathered in various different mass bins but coming from 
the same experiment. Given a set of effective coefficients g for all mass bins considered, one can 
investigate its dependence on Q. There is evidence that g increases slowly with Q, and in fact 
it is possible to fit a two-parameter function for the effective parameter g as given by (9), i.e. 
g = g 2 ln(Q/2Q )+g 1 , to this data set. We obtain g x = 0.15±0.13 GeV 2 , g 2 = 0.37±0.14 GeV 2 , 
with x 2 /d.o.f. = 0.33. These values are within the error range of the BLLY values (13) and lead 
to effective coefficients g for BLLY and our values of gi, g 2 being close to each other over the 
large range of Q, see Fig. 1. 

Due to the method of fitting, our errors could be larger then the ones listed in [25]. The 
normalization factors for the theoretical predictions with respect to the experimental data for 
the E605 and R209 experiments are found to be 0.91 ± 0.07 and 1.06 ± 0.09 respectively, and 
are in good agreement with [25]. The normalization of the CDF Z data from Run is kept 
equal to one, as in [25]. All our fitting is done using the MINUIT fitting programme [35] and 
MRST2001 parton distribution functions [36]. Given the differences in the fitting method and 
the parton distribution functions being used, we interpret our ability to recover (within errors) 
the previously determined BLLY parameters as an indication that our fitting method can be 
used as a valuable diagnostic tool. The small value of x 2 /d.o.f may suggest that the errors on 
the effective parameters g are overestimated. 

However, the BLLY fits were performed using data on Z production gathered by the CDF 
collaboration [31] during Run of the Tevatron pp collider. This particular data set has large 
statistical errors. The px distribution for Z production was measured much more accurately 
during Run I by both the CDF [33] and DO [34] collaborations. A new analysis of the form of 
the non-perturbative function, including the Run I Z data, appeared recently [26]. In agreement 
with the conclusions of [26], we find that the two-parameter form of the F NP ', given in Eq. 9, 
does not describe Run I data well, even after refitting for new values of parameters g\ and 
g 2 . In this analysis, in addition to CDF and DO Z data points with px < 20 GeV, we take 
R209 [30] data points in the 5 GeV < Q < 8 GeV bin only (p T < 2 GeV), E288 [32] data in 
the 5 GeV < Q < 6 GeV and 6 GeV < Q < 7 GeV bins (p T < 2 GeV, js = 27 A GeV) and 
E605 [29] data in the 7 GeV < Q < 8 GeV and 8 GeV < Q < 9 GeV bins also with p T < 2 
GeV. All points from CDF and DO Run I data sets have px < 20 GeV. Due to poor statistics, 
we do not include the Run data. Because of known normalization discrepancies between the 
CDF and DO experiments, normalization is allowed to be a free parameter while determining an 
effective parameter g separately for each data set. For data from other experiments, values of 
effective g parameters are determined in a fit which allows free normalization factor for theory 
predictions for each experiment, but these factors are kept the same for different Q bins from the 
same experiment. The normalization factors (multiplicating theoretical results) are: Ne 2 88 = 
0.81 ± 0.02, iV E6 05 = 0.90 ± 0.04, N R209 = 1.1 ± 0.1, N CDF = 1.09 ± 0.02, N D0 = 0.95 ± 0.03. 

x The inclusion of the 8 < Q < 11 GeV bin is dubious in any case, because of contamination by muon pairs 
coming from T decay. 
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Figure 1: Two-parameter fit to E605, R209 and CDF Z data (Run 0) data samples chosen as 
in [25] 

The fits of the b space predictions to data return a set of effective g coefficients for all mass bins 
considered, as shown in Fig. 2. A subsequent fit of the function g = g 2 In (Q/2Qo)+gi to this set 
returns g x = -0.08 ± 0.09 GeV 2 , g 2 = 0.67 ± 0.13 GeV 2 with x 2 /d.o.f = 8.25. Fig. 2 illustrates 
the difference between an effective coefficient g with these values of g±, g 2 and a coefficient g 
with the BLLY values (10) as the function of Q. 

The BLNY parameterization (15) assumes a gaussian form of the non-perturbative function 
which includes dependence on Q and x. To demonstrate this dependence in a clear fashion, it 
is convenient to rewrite the BLNY parametrization for the effective parameter g, cf. (15), 

9 = 91 + 92 In (t^) + 9193 ln ( 100x AXb) , (23) 
in the form 

9 = 9i+ » ln (4) +93hl (v / ?) ' <24> 

where we choose y/s = 19.4 GeV and Qo = 1.6 GeV. In effect, this parameterization introduces 
a y/s dependent component to g\ in the two-parameter form of the F NP in Eq. 9. Apart from 
the sets of data mentioned above we also include in the fits testing (24) the E288 [32] data in 
the 5 GeV < Q < 6 GeV and 6 GeV < Q < 7 GeV bins (p T < 2 GeV) for two additional cm. 
energies: y/s = 19.4, 23.8 GeV, as well as the E605 [29] data in the 10.5 GeV < Q < 11.5 GeV 
bin [pt < 2 GeV). All the data used in this analysis are contained in Table 1. Table 2 lists 
values of effective g coefficients returned by the fits to the data. They are also plotted in Fig. 3. 
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Figure 2: Two-parameter fit to E288, E605, R209, CDF and DO Z (Run I) data samples chosen 
as described in text. 

A subsequent fit of the function (24) to the set of these values reveals 

gi = 0.12 ± 0.07 GeV 2 , g 2 = 0.22 ± 0.12 GeV 2 , g 3 = 0.29 ± 0.09 GeV 2 , (25) 

with x 2 /d.o.f = 3.42. The best fit normalization factors (multiplicating theoretical results) are 
listed in Table 2. Obviously, the quality of the fit improves after introducing a third parameter. 
The quoted value of x 2 /d.o.f also suggests that the choice of logarithmic dependence on yfs in 
Eq. 24 could be too simplified. 

Translating the values of the g\, g 2 and g 3 parameters above to the BLNY non-perturbative 
function parameters (23) leads to 

g x = 0.26 ±0.1 GeV 2 , g 2 = 0.51 ±0.07 GeV 2 , g 3 = -0.55 ±0.011 GeV 2 , (26) 

These values (with the exception of g 2 ) lie within the error band of the parameters determined 
by the BLNY collaboration. Note, however, that we use a different fitting technique as well as 
different data samples for the analysis. 

We next turn our attention to the non-perturbative function for the case of Higgs production 
via gluon-gluon fusion, where the transverse momentum of the Higgs particle at small values of 
Pt is a result of soft gluon emission off the initial-state gluon lines. In our analysis we use E605 
data for T (pN) hadroproduction, allowing the normalization to be a free parameter fitted to 
data together with the non-perturbative parameters. We find the value of the effective parameter 
g(T) = 0.68 ±0.03 GeV 2 with x 2 /d.o.f = 2.79. Interestingly, this is relatively close to the values 
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Experiment 


(GeV) 


Q range (GeV) 


p T range (GeV) 


E605 


38.8 


7-9, 10.5-11.5 


0-2 


E288 


19.4, 23.8, 27.4 


5-7 


0-2 


R209 


62 


5-8 


0-2 


CDF 


1800 


91.2 


0-20 


DO 


1800 


91.2 


0-20 



Table 1: List of data used in the three-parameter analysis. 



Experiment 


(GeV) 


Q range (GeV) 


effective g (GeV^) 


Normalization 




E605 


38.8 


7-8 


0.414 ± 0.047 


0.842 ± 0.029 


3.126 






8-9 


0.472 ± 0.029 










10.5-11.5 


0.517 ± 0.038 






E288 


19.4 


5-6 


0.273 ± 0.028 


0.807 ± 0.020 


1.680 






6-7 


0.296 ± 0.031 








23.8 


5-6 


0.303 ± 0.029 










6-7 


0.357 ± 0.038 








27.4 


5-6 


0.354 ± 0.024 










6-7 


0.379 ± 0.146 






R209 


62.0 


5-8 


0.496 ± 0.071 


1.103 ± 0.096 


1.515 


CDF 


1800 


91.2 


2.746 ± 0.212 


1.090 ± 0.022 


0.509 


DO 


1800 


91.2 


2.615 ± 0.193 


0.954 ± 0.025 


0.911 



Table 2: Effective g coefficients from fits of b space theoretical predictions to data. 



of g for the corresponding qq Drell-Yan production process with invariant masses of the same 
order of magnitude, cf. Fig. 3. This conclusion is not unexpected, given the similar shapes of 
the pt distributions for Drell-Yan and T production data, see Fig. 4. 

The small difference between the values of g for Upsilon and Drell-Yan production suggests 
that the common prescription of multiplying the D-Y coefficient #2 by a factor Ca/Cf for gluon- 
initiated processes might require adjustments in the values of the other coefficients in order to 
maintain the quality of the fit. This is illustrated in Fig. 5, where apart from the effective values 
of g plotted for the analysed mass bins, we also plot the fitted effective g at \fs = \/^E605 = 38.8 
GeV, (23) and (24), with their coefficients #2, as listed in (16) and (25), multiplied by the 
factor Ca/Cf- The significant difference between the two predictions originates not only in 
the different values of the 52 coefficient obtained in our and the BLNY fits, but also from the 
different role of the gi coefficient itself - in our case (24) it solely measures the dependence on 
Q, whereas in (15) there is additional Q dependence in the coefficients g\ and gz through the 
relation x\%2 = Q 2 /s. As mentioned earlier, to bridge the gap between experiment and theory 
some other modifications to the values of the non-perturbative parameters will be needed in 
addition to the Ca/Cf prescription. In particular, one may argue that while the 53 parameter 
should stay the same in the non-perturbative functions relevant to processes with qq and gg 
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Figure 3: Three-parameter fit to E288, E605, R209, CDF and DO Z (Run I) data samples chosen 
as described in text. The parallel lines correspond to non-perturbative function of the form (24) 
with coefficients (25) and values of y/s of each experiment analysed. The line marked 'BLNY' 
corresponds to the BLLY fit of the form (15) with coefficients (16) at y/s = 38.8 GeV. 

initial states, the g\ coefficient can vary depending on the initial state. The universality of the 
53 coefficient can be motivated by general arguments for the logarithmic dependence on y/s for 
the contributions to the px of the final state from fluctuations in the underlying event. Figure 5 
shows the non-perturbative function (24) at y/s = \/*E605 = 38.8 GeV with g\ modified to 
obtain g = g(T) both when gi is multiplied by Ca/Cf and when it is left unchanged. The 
corresponding values of the 'adjusted' g\ parameter are g\ = —0.05 ± 0.4 GeV 2 for the former 
and g\ = 0.24 ±0.23 GeV 2 for the latter model. Consequently, the range of allowed values of the 
effective parameter g is rather large in each non-perturbative model discussed above. Moreover, 
since at the present time it is difficult to discriminate between the various models experimentally, 
the values of g from all reasonable models are allowed. A band of allowed values of the effective 
parameters g obtained in this way can be viewed as a useful estimate of the theoretical error in 
the pt distribution due to the non-perturbative input. 

4. Non-perturbative function in px space 

In general, the analysis of the non-perturbative function for the px space resummation method 
follows the b space analysis described in the previous section, except for the additional depen- 
dence of the non-perturbative function on the parameter pxiim- The value of this px space 
parameter is a priori unknown, and its determination requires simultaneous fitting alongside 
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Figure 4: E605 Drell-Yan (7 GeV< Q <8 GeV, 8 GeV< Q <9 GeV, 10.5 GeV< Q <11.5 
GeV bins) and T data compared to theoretical predictions using the best three-parameter fit. 
In order to better compare the shapes of the pt distribution, the Drell Yan data have been 
rescaled by the factor of 0.3, 0.6, 2.1 for the 7 GeV< Q <8 GeV, 8 GeV< Q <9 GeV and 10.5 
GeV< Q <11.5 GeV bins in Q, respectively. 



the other non-perturbative parameters. In b space, the analogous parameter bn m is traditionally 
fixed at a value of 0.5 GeV -1 . The dependence of the resummed distributions for Drell-Yan pro- 
duction on the parameter by im has been studied in [17]. In general, varying the b y m parameter 
will result in larger errors on the effective parameters g, and consequently on the transverse mo- 
mentum distribution itself. In the pt space case, all the fits we perform are fits for the effective 
a parameters (assumed different for each data set), the normalization factors (assumed identi- 
cal for data which come from the same experiment but have different Q) and the factor pxwm 
(assumed the same for all data sets). In fits for the pt space non-perturbative parameterization 
we use the same data sets and number of data points as in the b space fits. 

From a simultaneous fit to low Q Drell-Yan and Run I Z production data, we establish 
the value of priim to be PThm = 5.5 ± 2.98 GeV with x 2 /d.o.f = 1.54. The value of PThm is 
strongly correlated with the normalization factors for the low Q experiments and the value of 
a for Z production data, see Fig. 6. The normalization factors listed in Table 3 are different 
from their b space counterparts. The effective parameters a for all Q bins are now determined 
in one fit which also determines the common factor priim- In the b space analysis, the effective 
parameters g for each Q bin resulted from separate fits for all experiments - we did not fit for a 
common factor b\i m . The resulting value of PT\im appears rather large, if interpreted as defining 
the border between perturbative and non-perturbative physics. 
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Figure 5: Drell-Yan (E288, E605, R209), Z production (CDF and DO Run 1) and E605 Upsilon 
data together with the three-parameter fits, (15) and (24) at yfs = 38.8 GeV, with coeffi- 
cients (16) and (25) modified as discussed in text. 




Figure 6: The 68% C.L. contour in the priiim & cdf plane. 



Motivated by the form of the three-parameter non-perturbative function (24) in b space, we 
propose the following form for the dependence of the effective non-perturbative parameter a on 
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Q in (19): 

H^K^H'-G/sr <27) 

As for the b space analysis, we first fit the theoretical predictions to data in order to determine 
the effective coefficients a for all mass bins analysed. A fit of (27) to the set of effective values 
of a, listed in Table 3 and shown in Fig. 7, returns the following values for the parameters: 

fii = 0.20 ±0.50 GeV~ 2 , a 2 = 0.95 ± 0.92 GeV -2 , a 3 = 1.56 ± 0.57 GeV~ 2 , 

(28) 

with x 2 /d.o.f = 1.38. The errors on the coefficients are substantially larger than the correspond- 
ing errors on the coefficients in b space, again due to the much larger number of parameters (i.e., 
the effective values of a for bins in Q) determined in a single fit in the pt space analysis as 
compared to the b space fits. In Fig. 7 we also plot the non-perturbative parameterization (27) 
for the of each experiment analysed here, and the best fit coefficients oi, 02, 03 listed above. 
As for b space, the three parameter fit provides a much better description of the data than the 
two parameter case. 



Experiment 


(GeV) 


Q range (GeV) 


effective a (GeV 2 ) 


Normalization 


X 2 /d.o.f 


E605 


38.8 


7- 8 

8- 9 
10.5-11.5 


0.645 ± 0.149 
0.558 ± 0.072 
0.523 ± 0.079 


0.928 ± 0.067 


1.552 


E288 


19.4 


5- 6 

6- 7 


1.182 ± 0.216 
1.014 ± 0.219 


0.991 ± 0.050 






23.8 


5- 6 

6- 7 


0.998 ± 0.149 
0.865 ± 0.149 








27.4 


5- 6 

6- 7 


0.773 ± 0.080 
0.721 ± 0.047 






R209 


62.0 


5-8 


0.529 ± 0.252 


1.141 ± 0.288 




CDF 


1800 


91.2 


0.070 ± 0.012 


1.050 ± 0.065 




DO 


1800 


91.2 


0.070 ± 0.013 


0.914 ± 0.072 





Table 3: Effective coefficients a from fits of the pr space theoretical predictions to data. 

Given the b space results for the fit to the T data, we would expect the effective parameter a 
for T production to be relatively close to the corresponding values of a for Drell-Yan production. 
This is indeed the case, as can be seen from Fig. 7. The fitted a value for Upsilon production is 
o(T) = 0.39 ± 0.02 GeV~ 2 {x 2 /d.o.f. = 2.70). The value of o(T) is determined from a simple 
fit to the pt distribution data where the only other quantity fitted is the overall normaliza- 
tion; pTUm is kept fixed at the value determined from the fits to the DY data and the error 
on pt iim is not taken into account. The resulting error on the fitted value of a(T) is therefore 
superficially small. Nevertheless, it remains true that the central value of the effective a(T) is 
closer to the value predicted by the non-perturbative parametrization (27) with (28) than by 
the same parameterization with the 02 coefficient multiplied by Ca/Cf - the equivalent of the 
b space prescription of rescaling the g 2 coefficient by Ca/Cf- It therefore becomes an attrac- 
tive alternative to consider other possible modifications to the values of the non-perturbative 
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Figure 7: Effective parameters a for the Drell-Yan (E288, E605, R209), Z production (CDF and 
DO Run 1) and E605 Upsilon data and the non-perturbative parameterization (27), together 
with the best fit values of coefficients (28) plotted for the y/s of each set of experimental data 
analysed. 



parameters, determined for the Drell-Yan type processes, to describe processes with gluons in 
the initial state. Following the b space analysis presented in the previous section, we make the 
choice of adjusting the at parameter while keeping unchanged and the parameter a,2 either 
unchanged or rescaled by Ca/Cf- The resulting values of the parameter 2i are: 5i = 0.47l}'52 
GeV -2 and a\ = —0.82^1'^ for the former and latter model, respectively. The magnitude of 
the errors is a straightforward consequence of the large errors on the effective parameters a for 
the bins in Q. 



5. The Higgs transverse momentum distribution in the px space 
formalism 

As discussed previously, in the resummation framework the pr distribution for the gluon-fusion 
induced Higgs production process is obtained assuming that the non-perturbative contribution 
has the same form as the non-perturbative function for Drell-Yan type processes — but with 
the (72 or a corresponding parameter rescaled by a factor of Ca/Cf- The results presented in 
the previous two sections may suggest exercising some caution when applying this hypothesis. 

In this section we focus on the effect of the non-perturbative function on predictions for 
the Higgs pt distribution, in both the b space and pt space approach. In addition, we discuss 
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the changes made to the coefficients of the non-perturbative function(s) to allow for gluons 
in the initial state and the impact of these modifications on the predictions for the Higgs pr 
distribution. 

First, we focus on the b space non-perturbative function. As can be seen from Fig. 8, the 
central values of the effective parameters g for ^/s = 14 TeV are rather close to each other, 
irrespective of the different treatments of the gi, g2, g% coefficients discussed above. However, 
the significant errors on the coefficients g± , gi an d gs result in an even more pronounced error on 
the effective parameter g in each model. The errors are larger in models in which the 52 coefficient 
is rescaled by Ca/Cf, as the errors get rescaled together with the central values. If we assume 
the mass of the Higgs boson to be Mh = 125 GeV, the effective values are: g = 3.82 ±1.66 GeV 2 
for the model with 52 rescaled by Ca/Cf and g±, g% unchanged, g = 3.65 ± 1.99 GeV 2 for the 
model with 52 rescaled by Ca/Cf, <?i modified and g% unchanged, g = 2.96 ± 1.25 GeV 2 for 
the model with (72 , 93 unchanged and g\ modified. The next step is to implement the non- 
perturbative function (22) with these effective values of g into the resummed formalism. Again, 
we view the spread in the central values of g due to various models, together with the error 
bands on g in each model, as a measure of the theoretical error on the effective parameter g. 
Consequently, we present the Higgs px distributions for the smallest and the largest value of g 
allowed by the error bands for all models. These values both happen to come from the model 
with #2 rescaled by Ca/Cf, <?i modified and g% unchanged, and correspond to g = 1.67 GeV 2 
and g = 5.64 GeV 2 . 

Figure 10 shows the predicted transverse momentum distribution of the Higgs boson at the 
LHC calculated in the framework of the standard b space resummation with the effective non- 
perturbative parametrization of the form (22) for two 'extreme' values of g listed above. The 
resummed predictions are accurate up to NNLL level, i.e. they include the NNLL coefficients 
CW, as W ell as a numerical estimate of the coefficient [37], rescaled by Ca/Cf- 

The matching is performed at the leading-order level. We assume Mh = 125 GeV and use 
MRST2001 parton distribution functions. We find, in agreement with Ref. [10, 11], only a 
very small dependence of the theoretical predictions on the exact values of the non-perturbative 
parameters. The predictions shown in Fig. 10 are somewhat lower than the ones presented in [12] 
but higher then those obtained in [11]. 

Analogously to b space, the central values of the p^-space parameter a, as predicted by 
the different models for ^/s = 14 TeV and Mh = 125 GeV, lie relatively close to each other, 
see Fig. 9. However, due to the large errors on the central values of a associated with each 
model, the range of possible values of a is an order of magnitude large. The central values (for 
M H = 125 GeV) are: a = 0.0541°;°^ GeV~ 2 for the model with a 2 rescaled by C A /C F and a 1 , 
03 unchanged, a = 0.058^QQ2g GeV~ 2 for the model with 02 rescaled by Ca/Cf, ai modified 
and 03 unchanged, a = 0. 070^026 GeV~ 2 for the model with 02, 03 unchanged and a\ modified. 
The smallest and largest values of a also happen to come from the model with a.2 rescaled by 
Ca/Cf, a\ modified and 03 unchanged, and are: a = 0.032 GeV~ 2 and a = 0.308 GeV" 2 . 

The Higgs transverse momentum distribution obtained in the framework of the pt space 
formalism is compared to the b space results in Fig. 10. The px space distribution is calculated 
assuming the effective non-perturbative function of the form (19) with a coefficient a = 0.06 
GeV -2 and the best fit value of px\im = 5.5 GeV. The value of a = 0.06 GeV~ 2 is an average 
of the central values of a in the various models of treating the d\ , 0,2, 03 coefficients discussed 
above. The px space prediction agrees reasonably well with the b space prediction in the small 
Pt regime. The gap between the locations of the peaks of the distributions gets narrower for 
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Figure 8: The effective parameter g as a function of Q in various models of the non-perturbative 
parametrization in b space at y/s = 14 TeV. The central values are denoted by thick lines; the 
errors bands by thin lines. The dotted vertical line is drawn for Q = 125 GeV. 

smaller values of pram as well as smaller values of a. At larger values of px, the predictions of 
the two formalisms differ substantially. It is well known that the b space distribution, even after 
matching with the fixed-order result, is eventually going to turn negative at sufficiently large pt- 
This happens because the matching relies on the perturbative fixed-order result, known up to a 
certain order in as- At one order higher than the matching accuracy, the resummed expression 
introduces logarithmic terms which will not be cancelled by the terms coming from the matching 
term Y in (6). Since here matching is performed only at the LO level, the negative behaviour of 
the matched resummed distribution at large pt is quite pronounced. In the numerical realisation 
of the pt space method, only a subset of subleading higher order terms is summed. Consequently, 
the matched resummed distribution becomes negative at much larger values of pr, increasing 
the range of applicability of the resummation approach. 

As the analytic results for the NLO pt Higgs distribution are available in the literature [38] , it 
is possible to improve the accuracy of results presented here by performing matching at the NLO. 
However, our main goal here is to (re-)examine the influence of the non-perturbative function 
on the pt distribution and for that purpose we consider the LO matching sufficient. The NNLL 
resummed b space resulted matched to the NLO fixed-order expression was recently obtained 
in [12]. The problem of the negative large pt behaviour is normally solved by simply switching 
from the matched resummed cross section to the fixed order cross section (at an appropriate 
value of pt) to obtain the cross section valid at all pr- 

At low values of pt, the b space predictions for the two choices of g values are close to 
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Figure 9: The effective parameter a as a function of Q in various models of the non-perturbative 
parametrization in px space at y/s = 14 TeV. The central values are denoted by thick lines; the 
errors bands by thin lines. The dotted vertical line is drawn for Q = 125 GeV. 

each other, in agreement with the conclusions of Ref. [10, 11]. The effect of changing the non- 
perturbative part is felt, albeit slightly, at all values of px- In contrast, all the predictions 
from the px space formalism are already equal at px ~ 15 GeV. Around this value, the px 
space formalism begins to return a 'pure' perturbative part (the Sudakov factor) without any 
non-perturbative contamination. It is clear from Eq. 17, that for px 3> PThm and in the limit 
a —* oo the px space expression for the resummed cross section has only purely perturbative 
content. In the region of px where the px space predictions differ, the dependence on the values 
of the effective non-perturbative parameter a is, however, much stronger than the analogous 
dependence on the parameter g in b space. For lower values of a than the central value a = 
0.06 GeV -2 , the position of the peak and the value of the distribution at the peak changes 
significantly. In the case of the lower value of a = 0.032 GeV -2 the value of the distribution at 
the peak is around 5% larger than the corresponding value for a = 0.06 GeV -2 and the peak 
location moves by about 4 GeV towards smaller values of px- For values of a larger than the 
central value, the difference in the shape of the distributions is not so significant at moderate 
Pt- 2 

We believe that the difference between the 6-space and py-space predictions for the central 
values of the non-perturbative coefficients in the small px region originates in the essential 

2 As an artefact of the method of implementing the non-perturbative contribution in the pr space formalism, 
the distribution becomes unphysical for large values of a (corresponding to a strongly peaked distribution) at very 
small pt- 
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differences in the methods for accounting for non-perturbative effects in the two formalisms. In 
particular, the b space formalism requires a prescription to regularize the argument of as at large 
b. This is achieved by introducing the quantity b* as in Eq. 7. The pr space formalism, in turn, 
does not require any regularization prescription. The error band on the predictions of the Higgs 
Pt distribution in the pt formalism, obtained by considering distributions for the two 'extreme' 
values of a, i.e. a = 0.032 GeV -2 and a = 0.308 GeV -2 , is much larger than the corresponding 
error band from the b space method. In particular, the lower value of a denotes the broadest 
non-perturbative gaussian component allowed, so that the contribution from the perturbative 
part is minimised in comparison with the other (allowed) values of a. However, it needs to be 
stressed that apart from intrinsic differences in the formalisms the methods of determining g 
and a coefficients for the Drell-Yan processes differ — to obtain the effective a coefficients we 
varied pT\\m which resulted in relatively large errors on the values of the determined parameters. 
It is possible that a more accurate determination of the parameters would reduce the size of the 
error band on the predictions for the Higgs pr distribution in the px space formalism. It is also 
true that allowing b\\ m to vary while determining the g parameters would lead to a larger error 
band on the predictions from the b space formalism. Nevertheless, the larger error band on the 
Pt space results may suggest that the corresponding theoretical uncertainty on the predictions 
obtained with the help of the b space formalism is underestimated. 

An alternative method for determining the values of the parameters for the px space non- 
perturbative ansatz (27) relies on fixing pT\\m while fitting data for the effective parameters 
a. This method of finding the non-perturbative coefficients for the pt space parameterization 
corresponds in a much more direct way to the b space analysis and results in a significant 
decrease in the errors bars on the value of the coefficients. However, when applied it returns a 
considerably larger value of x 2 /d.o.f for the fits for the a\, ai and parameters. 

6. Conclusions 

We have performed detailed fits of the resummed theoretical predictions to the Drell-Yan data 
on the pt distributions for both the low Q fixed target experiments and the Z boson distribution 
data from the Tevatron pp collider. From these fits we established values of the non-perturbative 
components present in the resummed expressions. The fits were performed for two resummation 
formalisms: the well-known b space formalism and the pt space formalism. In both cases 
the non-perturbative parametrization is assumed to take the form of a simple gaussian with a 
logarithmic dependence on Q and ^/s. 

The resummed expressions for the Higgs pt distribution in the gluon-gluon fusion production 
process also require a parametrization of non-perturbative effects. However, one may expect that 
the gg initial state relevant for Higgs production will have a different non-perturbative function 
supplementing the resummed expressions compared to the qq Drell-Yan process. We used a com- 
bination of theoretical intuition and a fit of the resummed distribution to data on T production 
in hadron-hadron collisions (which is also predominantly characterized by the gg initial state) 
to determine the values of the coefficients in the corresponding non-perturbative functions. To 
summarize, we see a difference in the peak region of the Higgs transverse momentum distribution 
between the b and pt space predictions. In the latter case, the non-perturbative dependence 
is significant at lower pt but dies off rather rapidly with pr- For the b space prediction, in 
contrast, the non-perturbative variation is smaller at low pt but more persistent at higher pt- 
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Figure 10: Higgs boson px distributions at the LHC, as predicted by the b space and px space 
resummation formalism. 



We believe that we are unable to rule out any of the predictions described in this study on the 
basis of theory and experiment, and therefore that the corresponding span in the predictions is a 
measure of the true uncertainty on the theoretical prediction for the small px Higgs distribution. 
We note also that the uncertainty in the Higgs px distribution due to the non-perturbative 
component in the px space method seems to be larger than in the b space method. 
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